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Stealthy potentials, a family of long-range isotropic pair potentials, produce inhnitely degen¬ 
erate disordered ground states at high densities and crystalline ground states at low densities in 
d-dimensional Euclidean space In the previous paper in this series, we numerically studied the 
entropically favored ground states in the canonical ensemble in the zero-temperature limit across 
the first three Euclidean space dimensions. In this paper, we investigate using both numerical and 
theoretical techniques metastable stacked-slider phases, which are part of the ground-state mani¬ 
fold of stealthy potentials at densities in which crystal ground states are favored entropically. Our 
numerical results enable us to devise analytical models of this phase in two, three and higher dimen¬ 
sions. Utilizing this model, we estimated the size of the feasible region in configuration space of the 
stacked-slider phase, finding it to be smaller than that of crystal structures in the inhnite-system-size 
limit, which is consistent with our recent previous work. In two dimensions, we also determine exact 
expressions for the pair correlation function and structure factor of the analytical model of stacked- 
slider phases and analyze the connectedness of the ground-state manifold of stealthy potentials in 
this density regime. We demonstrate that stacked-slider phases are distinguishable states of matter; 
they are nonperiodic, statistically anisotropic structures that possess long-range orientational order 
but have zero shear modulus. We outline some possible future avenues of research to elucidate our 
understanding of this unusual phase of matter. 

I. INTRODUCTION 

A fundamental problem of statistical mechanics is the 
determination of the phase diagram of interacting many- 
particle systems. A substantial variety of pair interac¬ 
tions can produce a dramatic diversity of macroscopic 
phases, including crystals [T], quasicrystals HHS], liquid 
crystals [7], hexatic phases [SHII], disordered hyperuni¬ 
form systems and liquids nn. While crystals 

and liquids are the most common condensed states of 
matter, there are other states in between. For example, 
quasicrystals and liquid crystals both have anisotropy 
and long-range orientational order, like crystals, but lack 
long-range translational order, similar to liquids. Other 
phases with features that lie between crystals and liquids 
include disordered hyperuniform systems, which are dis¬ 
ordered but behave more like crystals in the way in which 
they suppress long-range density fluctuations mui]. 

A family of long-range isotropic pair potentials, called 
stealthy potentials, produces infinitely degenerate disor¬ 
dered hyperuniform classical ground states at high den¬ 
sities in d-dimensional Euclidean space [HUHlIini- 
HSl US]- Stealthy potentials are often specially con¬ 
structed such that finding a ground state is equivalent 
to constraining the structure factor S'(k) to be zero for 
all wave vectors k such that 0 < |k| < K, where K 
is some radial cutoff value. A dimensionless measure of 
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the relative fraction of constrained degrees of freedom 
(proportional to AT^) compared to the total number of 
degrees of freedom, x, controls the degree of order and 
degeneracy of the ground states of these potentials. 

In the preceding paper [^, we numerically studied 
the entropically favored ground states, i.e., configurations 
most likely to appear in the canonical ensemble in the 
zero-temperature limit, of stealthy potentials. We found 
that entropically favored ground states are disordered for 
X < 1/2, and crystalline for x > 1/2 up to a certain 
critical value m- 

The main focus of this paper is the investigation of 
stacked-slider phases, which are metastable states that 
are part of the ground-state manifold for some x above 
1/2, although not entropically favored. Stacked-slider 
phases were first discovered in two dimensions in Ref. [20] 
and were originally called wavy crystals because they 
were observed to consist of particle columns that display 
a meandering displacement away from linearity. How¬ 
ever, we will see that “stacked-slider phases” for arbi¬ 
trary dimensions is a more suitable name for this phase 
and this designation will be used henceforth. 

The authors of Ref. [501 easily distinguished stacked- 
slider phases from crystal phases by a lack of periodicity 
in direct space and a lack of Bragg peaks in its diffrac¬ 
tion pattern. Distinguishing stacked-slider phases and 
disordered phases, on the other hand, was based on a 
different property. In disordered phases, all k’s such that 
|k| > AT have positive structure factors. However, in 
stacked-slider phases, the structure factor at some k’s 
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such that |k| > K are implicitly constrained to vanish 
identically [20], i.e., they are induced to be zero by the 
constraints inside the radius K. The existence of implicit 
constraints was used to distinguish stacked-slider phases 
from disordered phases in Ref. [201 

There are still many outstanding questions concern¬ 
ing stacked-slider phases. Can a theoretical model of 
stacked phases in the thermodynamic limit be devised 
to elucidate previous numerical studies? One disadvan¬ 
tage of numerical studies is that finite-size effects make it 
difficult to conclude anything definitive about the large 
system limit. For example, are stacked-slider phases 
isotropic or anisotropic in this limit? Moreover, to what 
extent does the choice of the simulation box shape af¬ 
fect the results? Were any important features of stacked- 
slider phases overlooked by studying finite-precision sim¬ 
ulation results? Finally, because Ref. studied two di¬ 
mensions only, we do not know whether stacked-slider 
phases exist in other dimensions. This paper provides 
additional insights into these unanswered questions. 

The rest of the paper is organized as follows. In Sec. [IT] 
we perform numerical studies with much higher preci¬ 
sion than previously. The numerical results enabled us 
to find an analytical model of two-dimensional stacked- 
slider phases, presented in Sec. m We generalize this 
model to higher dimensions in Sec. |IV| We demonstrate 
that stacked-slider phases are distinguishable states of 
matter; they are nonperiodic, statistically anisotropic 
structures that possess long-range orientational order but 
have zero shear modulus. The model also shows that im¬ 
plicit constraints exist. In Sec. |V|we use this analytical 
model to show that stacked-slider phases are not entropi- 
cally favored in the zero-temperature limit of the canoni¬ 
cal ensemble. In Sec. |VI| we postulate that the transition 
between stacked-slider phases and disordered phases oc¬ 
curs at a slightly lower y than that reported in Ref. [20] 
from energy minimizations from high-temperature limit 
(Poisson) initial configurations. In Sec. |VII| we make con¬ 
cluding remarks and draw comparisons to other common 
phases of matter. 


II. NUMERICAL STUDY OF 2D 
STACKED-SLIDER PHASES 


In this section we numerically study the ground states 
of a stealthy potential at a variety of y’s (or densi¬ 
ties) in two dimensions. We begin with the mathemati¬ 
cal relations and simulation procedure in Sec. IIA and 
then present our results in Sec. El These results will 
suggest an analytical model of two-dimensional stacked- 
slider phases in Sec. |III[ 


A. Mathematical relations and simulation 
procedure 


As detailed in the preceding paper |2S| and other ref¬ 
erences [H nannHia, we simulate systems consisting 
of N point particles, located at = ri, r 2 , ..., r^r, in 
a simulation box in under periodic boundary condi¬ 
tions. The number density is p = N/vp, where vp is 
the volume of the simulation box. The particles inter¬ 
act with a pairwise additive potential u(r) such that its 
Fourier transform is: 


u(k) 


U(|k|) if |k| < K 
0 otherwise. 


( 1 ) 


where 'O(k) = f v(r) exp(—zk • r)dr is the Fourier trans¬ 
form of the pair potential u(r), V{k) is a positive func¬ 
tion, and AT is a constant. 

Under such potential, the total potential energy of the 
system can be calculated in the Fourier space 

= + ( 2 ) 

k 

where the sum is over all reciprocal lattice vector k’s 
of the simulation box such that 0 < |k| < K, n(k) = 
exp(-zk • Tj), and 


<^>o = [N{N-1) - N'^i(k)]/2vp (3) 

k 

is a constant independent of the particle positions r-^. 
Thus, the first term on the right-hand side of Eq. ([^ 
is the only configuration-dependent contribution to the 
potential energy 

= ^ E ^(k)l^^(k)r (4) 

^ 0<|k|<i<' 


Since V{k) >0 and vp > 0, Eq. Q shows that $*(r^) > 
0. Therefore, if configurations such that $*(r^) = 0 ex¬ 
ist, then they are the classical ground states of this po¬ 
tential. These configurations are achieved by constrain¬ 
ing n(k) to zero for all 0 < |k| < K and are said to be 
stealthy up to K. Since n(k) is related to the structure 
factor S'(k) by S'(k) = \h{'k)\‘^/N for every k 7^ 0, con¬ 
straining n(k) to zero is equivalent to constraining iS'(k) 
to zero. Let M be half the number of k points in the 
summation of Eq. ([^ [24| ; the parameter 


M 

d{N -1) 


( 5 ) 


determines the degree to which the ground states are 
constrained and therefore the degeneracy and disorder 
of the ground states [20|. For a fixed AT, the parameter 
X is inversely proportional to the density na mi US). 
When X ^ Xmaxi where Xmax is a dimension-dependent 
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constant, all the constraints are indeed satisfiable, thus 
<i)*(r^) of the classical ground states is zero [18]. The x 
values we study in this paper are always less than xtnax- 

In this section we choose N = 100, K = 1, and 
V{k) = 1. The relatively small choice of N increases 
the precision of the ground states we find. We will see 
that high precision is important in extracting an analyt¬ 
ical model from numerical results. The constant K and 
the magnitude of V(k) simply set the length scale and the 
energy scale. Although the function form of V{k) could 
theoretically affect the probability of sampling different 
parts of the ground-state manifold, it does not affect the 
manifold itself naiis]. As explained in Ref. [23 we use 
a rhombic simulation box with a 60° interior angle to 
alleviate finite-size effect. 

The ground states reported in this section are produced 
by the following steps. 

1. Start from a Poisson (i.e., ideal gas) initial config¬ 
uration. 

2. Minimize $*(r^) [in Eq. @] using the low-storage 
BEGS algorithm [25H25] . 

3. Minimize <I>*(r^) using the MINOP algorithm [25] . 

4. If $*(r^) < 10“^°, we successfully find a relatively 
high-precision ground state. 

5. Otherwise, what we find is either an imprecise 
ground state or a local minimum of $*(r^). There¬ 
fore, we discard this configuration. 

As detailed in Ref. [23 the low-storage BEGS algorithm 
is the fastest in minimizing $*(r^), while the MINOP al¬ 
gorithm finds the most precise ground states. Therefore, 
we minimize $*(r^) using these two algorithms consec¬ 
utively to maximize both efficiency and precision. 

These steps are performed Nt times for a variety of 
simulation box side lengths (and therefore a variety of 
X’s), listed in Table [l] [SD]- As detailed in Ref. (25] for 
a finite system, only certain values of x allowed. 
The X values in Table [^ contain all possible choices in 
the range 0.5 < x < xtnax^ which covers the previously 
reported stacked-slider phase regime in two dimensions, 
0.57... < X < 0.77... [23122]. Except for x = 0.8787..., 
where we could not precisely identify ground states, we 
plot the real-space configuration and reciprocal-space 
structure factor of at least 50 successful energy minimized 
results and visually inspect them. We divide them into 
different categories based on their appearances and then 
present representative configurations below. 

B. Results 

Representative numerically obtained ground-state con¬ 
figurations and their structure factors (in logarithmic 
scales) are presented in Eigs. [I]|^ For 0.5303... < x ^ 
0.6060 ..., the ground-state manifold appears to contain 
a variety of structures (see Fig. [^. Except for the first 
one, all real-space configurations in Fig. [^ appear to be 


TABLE I. The x values, number of trials Nt, and number of 
successes Ns for each simulation box side length L. 


L 

X 

Nt 

Ns 

56 

0.5303... 

1000000 

15615 

57 

0.5606... 

199915 

17127 

58 

0.5909... 

200000 

411 

59 

0.6060... 

199965 

8875 

60 

0.6363... 

1000000 

27788 

62 

0.6666... 

1000000 

76727 

63 

0.6818... 

1000000 

157501 

64 

0.7121... 

200000 

119563 

65 

0.7424... 

1000000 

165203 

66 

0.7575... 

200000 

80258 

68 

0.7878... 

200000 

2577 

70 

0.8787... 

200000 

0 


Bravais lattices. However, their structure factors are not 
as simple as a collection of Bragg peaks among a zero- 
intensity background, suggesting that the real-space con¬ 
figurations are not perfect Bravais lattices. 

At X = 0.6363... and x = 0 . 6666 ... a type 
of relatively-simple-looking configuration appears (see 
Fig.§. The real-space configurations appear to be com¬ 
prised of straight lines of particles with wavelike displace¬ 
ments relative to each other. The structure factors, on 
the other hand, consist of straight lines of nonzero values 
in a background of virtually zero (< 10 “^°) intensities. 

For X ^ 0.6818..., the results are similar to that 
in Fig. [^ but there exist so many constraints that the 
nonzero-value lines in the structure factor have to be 
interrupted. The interruptions grow in length as x 
increases and eventually, at x = 0.7878..., the only 
nonzero structure factors are the Bragg peaks and the 
real-space configuration becomes a Bravais lattice. 


III. ANALYTICAL MODEL OF 
TWO-DIMENSIONAL STACKED-SLIDER PHASE 

In this section we look closer at the simulation results 
that yield stacked-slider phases to see if an exact analyt¬ 
ical construction can be extracted. We will see that un¬ 
derstanding the configuration shown in Fig. [^is the key 
to understanding other configurations. The real-space 
configuration in Fig.[^seems to be made of straight hori¬ 
zontal lines that are displaced relative to each other. Are 
the displacements of different horizontal lines indepen¬ 
dent of each other or correlated in some way? To answer 
this question, we numerically constructed a configura¬ 
tion that is made of horizontal straight lines of particles, 
just like the one shown in Fig. [^ but with independent 
random displacements along each horizontal line. The 
structure factor of the new configuration has exactly the 
same support [the set of k’s such that S'(k) 7 ^ 0] as the 
one shown in Fig. [^ Thus, the new configuration is also 









4 



FIG. 1. (Color online) Four representative numerically obtained ground-state configurations at x = 0.5606... (left) and their 
corresponding structure factors (right), where colors indicate intensity values at reciprocal lattice points. 
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FIG. 2. (Color online) A numerically obtained ground-state configuration at x = 0.6363... (left) and the corresponding 
structure factor (right), where colors indicate intensity values at reciprocal lattice points. 
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FIG. 3. (Color online) Shown on the left are four representative numerically obtained gronnd-state configurations at x = 
0.6818 ... (first row), x = 0.7121.. . (second row), y = 0.7424 ... (third row), and x = 0.7878 ... (fourth row). On the right are 
their corresponding structure factors, where colors indicate intensity values at reciprocal lattice points. 


a ground state at this x value. Therefore, the displace¬ 
ments of each line do not need to be correlated in any 
way. This allows us to find a two-dimensional stacked- 
slider phase model, depicted in Fig. 


y 



FIG. 4. (Color online) Schematic plot of the two-dimensional 
stacked-slider phase model. Each horizontal line of particles 
[indicated by large (blue) dots] form a one-dimensional integer 
lattice with lattice spacing a. Then multiple horizontal integer 
lattices are stacked vertically, with spacing b. Each horizontal 
line of particles can be translated freely to slide with respect 
to each other. 


























































































































































































































6 


This analytical model allows the calculation of various 
properties of the two-dimensional stacked-slider phases. 
One can find the analytical pair correlation function and 
structure factor of this model, assuming that the dis¬ 
placement of each line is independent and uniformly dis¬ 
tributed between 0 and a. The pair correlation func¬ 
tion 32 (r) is defined such that pg 2 {v)dr is the conditional 
probability that a particle is found in the volume element 
dr about r, given that there is a particle at the origin. 
For the two-dimensional stacked-slider phase, 52 (r) can 
be found directly from the definition of this model: 

g 2 ix,y) = b'^5{y - bj) + ab5{y)'^5{x - aj), ( 6 ) 

3^0 

where x and y are horizontal and vertical coordinates, 
both summations are over all nonzero integers j, and <5 
denotes the Dirac delta function. The structure factor 
S'(k) can be found by Fourier transforming 52 (r) — 1: 

5(k) = l + pF[ 52 (r)-l], (7) 

where F[- ■ ■ ] denotes Fourier transform. Substituting ®) 
into Q, one gets: 

_ 27r5(fc„) l^27r ^ 

S[k.) — I , IIl27r/b(^y) 1 I 

a \o / ( 8 ) 

Ott 47r^ 

+-- T-5{kx)5{ky), 

where and ky denote the horizontal and vertical com¬ 
ponents of k, respectively, and IIIt (0 = Sj=^oo ~ 
jT) is the Dirac comb function. Both the pair correlation 
function ^ and the structure factor Q are anisotropic, 
since swapping x and y in (§ and ([^ gives different ex¬ 
pressions. 

A topological property this model can predict is the 
connectedness of the ground-state manifold, i.e., whether 
or not a ground state can be continuously deformed to 
another ground state without crossing any energy bar¬ 
rier. Each stacked-slider configuration is obviously con¬ 
tinuously connected to a rectangular lattice by the slid¬ 
ing motion of different lines. However, there are many 
permutations of the rectangular lattice. Are these per¬ 
mutations connected to each other through vertical and 
horizontal sliding motions? In the Appendix, we show 
that for a finite-size rectangular lattice consisting of N 
particles, all permutations are connected if and only if N 
is even. 

Having found an analytical model of the ground states 
in this X range, we move on to lower and higher x 
ranges. The lower x simulation results appear to be 
more complex. The first configuration in Fig. [l] ap¬ 
pears to be similar to our existing analytical model, 
except that the nonzero-value regions in the structure 
factor are not strictly lines: The highest-intensity lines 
[S'(k) ~ 10°] are surrounded by lower-intensity regions 
[S'(k) ~ 10“which are surrounded by even lower- 
intensity regions [S'(k) ~ 10“^°]. The structure factor 


in the lower-intensity regions are very small, but are 
still much larger than the machine precision. (We use 
double-precision numbers, which have around 16 signif¬ 
icant digits, to calculate h(k). Therefore, the machine 
precision of S{k) = \n(k)\'^/N should be on the order of 
(10“^®)^/Af = 10“^'^.) So a natural question arises: Are 
the lower-intensity regions real or are they an artifact of 
finite-precision simulations? 

To answer this question, we chose a k point right next 
to the highest-intensity line and plotted the structure fac¬ 
tor at this k point versus the potential energy during the 
energy minimization (see Fig. [^. As goes to zero, the 
structure factor at this k point also goes to zero. Thus, 
we believe the lower-intensity regions are the result of 
numerical imprecision. If one could carry out an infinite- 
precision simulation and drive this configuration to a true 
ground state, the structure factors in the lower-intensity 
regions should go to zero and the configuration would 
become consistent with our analytical model. 



FIG. 5. (Color online) (a) A numerically obtained ground 
state at y = 0.5606.... (b) The corresponding structure fac¬ 
tor. A specific k point is indicated by a black square and 
an arrow, (c) The structure factor at this particular k point 
is plotted against total energy 4>* during the optimization, 
showing S'(k) —>■ 0 as 'F* —>• 0. 

Having understood the first configuration in Fig. let 
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us move on to other configurations in that hgure. The 
second and third conhgurations appear to be intermedi¬ 
ate configurations between the first one and the fourth 
one. The fourth configuration looks like a Bravais lattice, 
except that the Bragg peaks are smeared out. Again, to 
find out whether this broadening of the Bragg peaks is 
real or artificial, we plotted the structure factor at a k 
point near a Bragg peak versus the potential energy in 
Fig.§ We find again that the structure factor at this k 
point goes to zero as goes to zero. Thus, the smearing 
out of the Bragg peaks is also due to numerical impreci¬ 
sion. If one could carry out an infinite-precision energy 
minimization on this configuration, one should get a Bra¬ 
vais lattice. 
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FIG. 6 . (Color online) (a) A numerically obtained ground 
state at y = 0.5606.... (b) The corresponding structure fac¬ 
tor. A specific k point is indicated by a black square and 
an arrow, (c) The structure factor at this particular k point 
is plotted against total energy il?* during the optimization, 
showing ^(k) —>■ 0 as >!>*—>• 0 . 


So far we have demonstrated that the numerically 
obtained ground states follow a simple model at y = 

0.6363 ... and y = 0.6666_We have also demonstrated 

that while the numerically obtained ground states for 
0.5303 ... < y < 0.6363 ... appear to be richer, they are 


actually exactly the same as either the model or a Bravais 
lattice if we could perform infinite-precision simulations. 
However, as we move to higher y’s, the ground states 
start to lose degrees of freedom. As shown in Fig. 
at y = 0.6818..., the high-intensity lines in the struc¬ 
ture factor develop zero-intensity interruptions. In our 
stacked-slider phase model, if each line of particles could 
move independently, then the high-intensity lines in the 
structure factor would have no interruptions. Thus, these 
interruptions indicate constraints in the displacements of 
each line of particles. At y = 0.7121..., the lines are in¬ 
terrupted even further, indicating even more constraints 
in the displacements of each line. At y = 0.7424..., the 
structure becomes a two-particle-basis crystal. Eventu¬ 
ally, at y = 0.7878..., the structure becomes a Bravais 
lattice. 

Starting from y = 0.6818..., the stacked-slider phase 
become more constrained as y increases. To study how 
constrained this phase is at different y values, we cal¬ 
culate the number of zero eigenvalues Ue of the Hessian 
matrix of the potential energy. This number is equal to 
the number of independent ways to deform the structure 
such that the energy scales more slowly than quadratic, 
which is an upper bound of the dimensionality of the 
ground-state configuration space Uc [i.e. the number of 
independent ways to deform the structure such that the 
$*(r^) remains zero]. For y < 0.6818 ..., our model pre¬ 
dicts ric = 11 (since there are two translational degrees 
of freedom, and nine independent ways to slide the ten 
lines of particles relative to each other) and our calcu¬ 
lation also find ng = 11. At y = 0.6818..., 0.7121..., 
0.7424 ..., and 0.7878 ..., our calculations find rig = 9, 5, 
3, and 2, respectively. This calculation suggests that as 
y increases, rig gradually decreases. Eventually, rig = 2, 
indicating that there is no way to deform the structure 
other than trivial translations. 


IV. GENERALIZED STACKED-SLIDER PHASE 
MODEL 

We now generalize the two-dimensional stacked-slider 
phase model to higher dimensions. To begin with, we 
present and prove the following theorem: 

Stealthy Stacking Theorem. Let dp and dg be posi¬ 
tive integers. Let W be (dp-|-dQ)-dimensional Euclidean 
space. Let Wp he a dp-dimensional subspace of W and 
Wq be the dg-dimensional orthogonal complement space 
of Wp. Let P be a point pattern in Wp with density pp. 
Eor each point a £ P, let Q{a) be a point pattern in Wq 
with some density pg independent of a. If P is stealthy 
up to certain reciprocal-space cutoff Kp and all (5(a)’s 
are stealthy up to certain reciprocal-space cutoff Kg in 
their subspace, then the following point pattern in IT, 

{a-f bja G P, b e Q(a)} (9) 

is a stealthy point pattern up to AT = min{Kp, Kg). 















Proof. The collective density variable of the point pat¬ 
tern in Eq. is 

= exp[-ik- (a-f b)]. (10) 

aGPbGQ(a) 

Since Wp and Wq are two orthogonal complementary 
subspaces of W, we can divide vector k into two parts 
k = kp-l-kg, where kp G Wp and kg G Wq. Therefore, 

^(k) = XI X exp[-i(kp-tkg) • (a-Gb)] 
aGPbGQ(a) 

= X®^P(*kp-a) ^ exp(-ikg-b). 

aGP bGQ(a) 

For any k such that 0 < |k| < K, |kg| < |k| < iF < Kq. 
If kg 7 ^ 0, then the stealthiness of point patterns Q(a) 
gives 


exp(-ikg • b) = 0 (12) 

bGQ(a) 

and therefore p.(k) = 0. On the other hand, if kg = 0, 
then kp = k and Eq. 0 becomes 


^(k) =iVQ(a) ^exp(-ik-a), (13) 

aGP 


where ./Vg^^) is the number of particles in pattern (3(a), 
which is independent of a because all the (3(a)’s have the 
same density. Since 0<|k|<iG<ifp, the stealthiness 
of point pattern P gives 

n(k) = 0. (14) 

To summarize, for any k such that 0 < |k| < if, whether 
or not kg = 0, h(k) is always zero. Therefore, the point 
pattern Q is stealthy up to K. 

The parameter x of this point pattern can be calcu¬ 
lated using Eq. (35) of Ref. [18]. Our calculation yields 

^ vi{dp + dQ;K) 

^ 2{dp + dQ){2Tr)‘^p+‘^Q ppPq' 


where vi{d;r) is the volume of a d-dimensional hyper¬ 
sphere of radius r. In the case Kp = Kq, using Eq. (35) 
of Ref. [T8|, Eq. (151 can be simplified to: 


X = 


2ui(dp-G dg; 1) dpdg 

ui(dp; l)ui(dg; 1) dp -G dg 


XPXQ- 


(16) 


The aforementioned theorem allows us to construct 
stacked-slider configurations in higher dimensions. To 
construct a stacked-slider configuration in d > 2, choose 
two lower dimensions dp and dg such that dp -G 
dg = d. Choose a dp-dimensional stealthy configu¬ 
ration P and replace each particle a in P with a dg- 
dimensional stealthy configuration (3(a) and the result¬ 
ing d-dimensional configuration is a stacked-slider one. 







FIG. 7. (Color online) Schematic plot of the stacked-slider 
phase model. The large black dots form an integer lattice 
(point pattern P). By replacing each black dot with a two- 
dimensional stealthy point pattern (indicated by small bine 
dots) of the same density [point patterns (3(a)], the over¬ 
all three-dimensional point pattern consisting of all the small 
blue dots is stealthy. The two vectors rp and rg are in sub¬ 
spaces Wp and Wq, respectively. Note that since some (3(a)’s 
are two-dimensional stacked-slider configurations, this config¬ 
uration allows both interlayer and intralayer sliding motions, 
as detailed in Sec. IIIII 


The resulting configuration is often anisotropic, since dp 
dimensions are treated separately from the remaining dg 
dimensions. See Fig. for an illustration of a three- 
dimensional stacked-slider configuration with dp = 1 and 

dg = 2. 

Certain three-dimensional crystal structures can allow 
sliding deformations while remaining stealthy at rela¬ 
tively large (greater than 0.5) %. As Fig. shows, the 
simple cubic lattice allows the sliding motion of each 
two-dimensional square-lattice layer and the sliding mo¬ 
tion of each line of particles inside every layer for x up 
to 0.6981.... Barlow packings |3T], including the face- 
centered-cubic packing and the hexagonal close packing, 
also allow the sliding motion of each triangular-lattice 
layer of particles for y up to 0.7600 .... 

Equation (161 can be used to calculate the maximum 
X values of the stacked-slider-phase, Xmaxi assuming un¬ 
constrained sliding motions, in each space dimension d. 
To do this one can try all possible combinations of pos¬ 
itive integers dp and dg such that dp -G dg = d, and 
let xp and XQ equal to Xmax in dp and dQ dimen¬ 
sions, respectively. Our calculations for 2 < d < 4 are 
summarized in Table El There is no obvious trend in 
these low dimensions. However, as d increases, the fac- 
toi- Eq. (@ decreases for any dp 

and dg. Thus, Xmax should become arbitrarily small in 
sufficiently high dimensions. 

Similar to two-dimensional stacked-slider configura¬ 
tions, the higher-dimensional ones also have implicit con¬ 
straints [i.e. k vectors such that |k| > AT and S'(k) = Oj. 
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FIG. 8. (Color online) Three-dimensional stacked-slider con¬ 
figuration stealthy up to x = 0.6981.... This configuration 
is obtained by sliding each vertical plane of particles relative 
to each other and then sliding each vertical line in each plane 
relative to each other starting from the simple cubic lattice. 


TABLE II. Comparison of the maximum y value of stacked- 
slider-phases predicted by the generalized model Xmax and 
the maximum x value of Bravais lattices Xmax in two, three, 
and four dimensions. 


d 

A.max 

\max 

\max /\max 

2 

tvIA 


0.8660... 

3 

471 

9\/3 

2\/27r 

9 

0.8712... 

4 

_ IS _ 

TT^ 

_8_ 

0.7071... 


As seen in Eq. (11), S'(k) = \n{'k)\^/N = 0 as long as 
0 < Ikgl < K. One can thus choose arbitrarily large kp 
such that |k| = |kp -|- kgl > K. 


V. FEASIBLE REGION OF THE 
CONFIGURATION SPACE 


states is much smaller than the entire feasible region in 
the configuration space, this set will almost never appear 
in the canonical ensemble, i.e., they are not entropically 
favored. 


In the infinite-system-size limit, the feasible region of 
any stacked-slider configuration is much smaller than 
that of any crystal if both the stacked-slider configura¬ 
tion and the crystal are ground states. This is because as 
N —>■ oo, the configurational dimension ric [i.e., the num¬ 
ber of independent ways to deform the structure such 
that the $*(r^) remains zero] of stacked-slider phases 
scales more slowly than the number of particles N. For 
example, for a two-dimensional stacked-slider configura¬ 
tion in which each row of particles can slide i nde oen- 


dently, ric scales as '/N- As discussed in Sec. Ill the 


number of zero eigenvalues of the Hessian matrix of the 
potential energy ng is equal to rig- Since a nonzero eigen¬ 
value of the Hessian matrix corresponds to a quadratic 
scaling in one direction, in the c?A^-dimensional configu¬ 
ration space, <I>*(r^) has quadratic scaling in dN — ric 
directions. In these directions, as e —> 0, the width of 
the feasible region scales as y/e. In the remaining rig di¬ 
rections, the width of the feasible region is much larger, 
since these directions correspond to translations of dif¬ 
ferent rows of particles, which keeps <i)*(r^) zero. If we 
let the widths of the feasible region in these directions 
be L, then the total volume of the feasible region of the 
stacked-sliding phase is approximately 


V « JJia^{dN-na)/2 


j^VN ^{dN-VN)/2 


(17) 


In the case of a crystalline structure, rie scales as N 
when N —i oo. This can be seen in Fig. where we plot 
/ = ric/dN versus N for triangular lattices at x = 0.6. 
This figure shows that / tends to some constant as N 
grows, which means rig scales as N. Since a zero eigen¬ 
value of the Hessian matrix of $*(r^) implies a slower- 
than-quadratic scaling in some direction, the width of the 
feasible region in these rig directions scales larger than ^/e 
as e —)■ 0. Let the widths of the feasible region in these 
rig directions be , where 0<a;<l/2is some exponent. 
The width of the feasible region in the remaining dN — ric 
directions scales as The total volume of the feasible 
region of a crystal is approximately 


Although stacked-slider configurations are part of the 
ground-state manifold of stealthy potentials, we will show 
in this section that they are not entropically favored, as 
indicated in Ref. |T51 Entropically favored ground states 
are the configurations that most likely appear in the 
canonical ensemble in the zero-temperature limit |18| . In 
this limit, as a good approximation, the system can only 
visit part of the configuration space where $*(r^) [in 
Eq. is less than e, where e > 0 tends to zero as the 
temperature tends to zero. This part of the configura¬ 
tion space is therefore called the feasible region. If the 
feasible region corresponding to one set of the ground 


V « ^(dN-nc)/2^ncX 


^dN(l-f)/2^dNfx 


(18) 


The ratio of 14 and 14 is approximately 
Yl ~ T\dN [dNf{l-2x)-VN]/2 


(19) 


Since x < 1/2, as A —>• oo and e —0, ^ 0. Therefore, 

the feasible region of the stacked-slider phase is much 
smaller than that of the crystal. Since there are always 
crystalline structures competing with the stacked-slider 
phase, the latter is never entropically favored. 
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FIG. 9. (Color online) Fraction of zero eigenvalues of the 

rie 

Hessian matrix of the potential energy / = for triangular 
lattices of various numbers of particles at x = 0.6. 


of the above discussion. Previously, the existence of im¬ 
plicit constraints [k’s such that |k| > K and <S'(k) = 0] 
was often used to distinguish stacked-slider phases from 
disordered phases [101 [H]. Therefore, one-dimensional 
stealthy ground states in the range 1/3 < x < 1/2, 
proven to have implicit constraints [32j . were considered 
to be stacked-slider phases [H]. However, this study sug¬ 
gests that one-dimensional stealthy ground states in this 
range are not a typical stacked-slider phase. First, our 
model only predicts stacked-slider phases if the space di¬ 
mension c? is a sum of two positive integers d = dp + dq. 
This requires that d > 2. Second, the x range of the 
one-dimensional stealthy ground states with implicit con¬ 
straints is also very different from that of the higher¬ 
dimensional stacked-slider phases. We also found that 
one-dimensional stealthy ground states in this x range 
satisfy the relation / = 1 — 2x and can be obtained from 
energy minimizations starting from random initial config¬ 
urations with 100% success rate; both are characteristics 
of disordered phases 


VI. RELATIVE STABILITY OF 
STACKED-SLIDER PHASES 


VII. CONCLUSIONS AND DISCUSSION 


We have shown that the feasible region of stacked- 
slider phases is always smaller than that of crystal phases 
and thus concluded that stacked-slider phases are never 
equilibrium phases at T = 0. This conclusion is con¬ 
firmed by low-temperature molecular dynamics simula¬ 
tions reported in Ref. [23 which found disordered struc¬ 
tures for X < 1/2 and crystalline structures for x > 1/2. 
However, this simple conclusion cannot explain or predict 
energy minimization results from high-temperature ini¬ 
tial configurations that were used previously [20] , where a 
transition from disordered phases to metastable stacked- 
slider phases was observed as x increases, characterized 
by the change of the support of 5'(k). In two dimensions, 
Ref.|20]reported that this transition is at x = 0.57..., but 
high-fidelity simulations, reported in Sec. jllj produced 
stacked-slider configurations at x = 0.5305..., suggest¬ 
ing that the transition is earlier than 0.5305 .... Another 
observation on the disordered region supports our result: 
Section V of Ref. reported that the fraction of normal 
modes with vanishing frequency / in disordered phases 
is exactly 1 — 2x for x < 1/2. However, this exact re¬ 
lation cannot be true for the x > 1/2 region, since / 
is non-negative. This suggests that there exists a sharp 
transition at x = 1/2, which is likely the phase transi¬ 
tion to the stacked-slider phase. Although Ref. [321 only 
reported the relation / = 1 — 2x in two dimensions, it 
explained this relation by simple counting arguments in¬ 
volving the number of constraints versus the number of 
degrees of freedom and hence this relation should apply in 
any dimension. Therefore, for any d, as long as stacked- 
slider phases exist for some x above 1/2, there should be 
a nonequilibrium phase transition from disordered phases 
to stacked-slider phases at the threshold x = 1/2. 

It is noteworthy that one dimension is an exception 


In this paper we studied using numerical and the¬ 
oretical techniques stacked-slider phases, which are 
metastable states that are part of the ground-state man¬ 
ifold of stealthy potentials at densities in which crystal 
ground states are favored entropically in the canonical 
ensemble in the zero-temperature limit [TH [21] . The nu¬ 
merical results suggested analytical models of this phase 
in two, three and higher dimensions. Utilizing this model, 
we estimated the size of the feasible region of the stacked- 
slider phase, finding it to be smaller than that of crystal 
structures in the infinite-system-size limit, which is con¬ 
sistent with our recent previous work [M123. In two 
dimensions, we also determined exact expressions for the 
pair correlation function and structure factor of the an¬ 
alytical model of stacked-slider phases, and analyzed the 
connectedness of the ground-state manifold of stealthy 
potentials in this density regime. 

Our analytical constructions demonstrate that 
stacked-slider phases are nonperiodic, statistically 
anisotropic structures that possess long-range orien¬ 
tational order but have zero shear modulus. Since 
stacked-slider phases are part of the ground-state mani¬ 
fold of stealthy potentials, they are also hyperuniform. 
Therefore, stacked-slider phases are distinguishable 
states of matter that are uniquely different from some 
common states of matter listed in Table [ml Note 
that distinctions between the attributes indicated 
in the table may be subtly different. For example, 
crystals, quasicrystals, and stacked-slider phases all 
have long-range orientational order, but with different 
symmetries. While crystals can only have twofold, 
threefold, fourfold, or sixfold rotational symmetries, 
quasicrystals have prohibited crystallographic rotational 
symmetries. Stacked-slider phases generally do not have 
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TABLE III. Comparison of the properties of some common states of matter. Here crystals and qnasicrystals signify perfect 
crystals and perfect qnasicrystals, respectively, without any defects (e.g., phonons and phasons). The checks and crosses indicate 
whether or not different phases have the attributes listed in the first column. 


Property 

Crystals 

[I] 

Quasicrystals 

[2H6] 

Stacked-slider 

phases 

Disordered 
ground states of 
stealthy potentials 

miiiiinHll! 

Liquid crystals 

m 

Liquids 

m 

periodicity 

/ 

X 

X 

X 

X 

X 

positive 
shear modulus 

/ 

/ 

X 

X 

X 

X 

hyp eruniformity 

/ 

/ 

/ 

/ 

X 

X 

anisotropy 

y 

/ 

y 

X 

/ 

X 

long-range 
orientational order 

/ 

/ 

y 

X 

/ 

X 


any rotational symmetry, but the fact that they can 
be constructed by stacking lower-dimensional stealthy 
configurations in a higher-dimensional space makes the 
stacking directions different from the sliding directions, 
giving them their unique orientational order. 

Our understanding of stacked-slider phases is only in 
its infancy with many open questions. For example, 
what is the nature of the associated excited states? Can 
stacked-slider phases emerge from particles interacting 
with other potentials not necessarily as ground states? 
Can such phases be entropically favored in some ensem¬ 
ble and with what other phases would it coexist? This is 
just a partial list of possible of future avenues of research 
in our understanding of this unusual phase of matter. 


am Ax B matrix. Three different sliding motions will be 
frequently used in this section. They are as follows: 


• Move the top row of particles to the right by one 
particle spacing, denoted by 

• move the leftmost column of particles upward by 
one particle spacing, denoted by =>; 


• and move the leftmost column of particles down¬ 
ward by one particle spacing, denoted by 

As an example of this notation, for A = i? = 2, per¬ 


mutations 



and 



are connected because 
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(Al) 


Similarly, permutations 
nected because 




are con- 


Appendix A: Connectedness of permutations of 2D 
stacked-slider phase 



(A2) 


As discussed in Sec. m each two-dimensional stacked- 
slider configuration is connected to a permutation of 
the rectangular lattice. Therefore, a natural question 
is whether or not these permutations of the rectangu¬ 
lar lattice are also connected through sliding motions. If 
all permutations of the rectangular lattice are connected, 
then the entire stacked-slider phase ground-state mani¬ 
fold is connected. We will show that, for a rectangular 
lattice consisting of A rows and B columns of particles, if 
each row and each column can slide individually, then all 
permutations of the rectangular lattice are connected if 
and only if AB is even. We will number all the particles 
from 1 to AB. Each permutation will be represented by 


So far we have demonstrated that it is possible to swap 
the two adjacent particles in the first row [by Eq. (All] 
or the two adjacent particles in the first column [by 
Eq. (A2|] for A = B = 2. Since the system has transla¬ 
tional symmetry, one can swap any two adjacent parti¬ 
cles. The swapping of any two nonadjacent particles can 
be done by a series of adjacent-particle swapping. For ex- 

/f 2^ 

ample, to swap nonadjacent particles 1 and 4 in I ^ 


one can swap particles 1 and 2, then swap particles 1 and 
4, and then swap particles 2 and 4. Finally, since we can 
swap any two particles, we can connect one permutation 
to any other permutation by swapping each particle with 
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the particle in its new place. Therefore, all permutations 
of 2 X 2 rectangular lattices are connected by row-sliding 
and column-sliding movements. 

Next, we show that one can swap two adjacent particles 
for A = 3 and 5 = 4. To swap the first two particles 
in the first row, one can perform the following sliding 
operations: 


A 

2 

3 

4\ 

(4 

1 

2 

3\ 


6 

7 

00 

5 

6 

7 

8 

V9 

10 

11 

12 J 

[9 

10 

11 

12/ 


/5 1 2 3\ 

9 6 7 8 

10 11 12 / 


/3 5 1 2\ A 5 1 2\ /2 4 5 l\ 

^ 96 7 8 ^ 36 7 8 ^ 36 7 8 

\4 10 11 12/ \9 10 11 12/ \9 10 11 12/ 


/3 4 5 l\ 

9 6 7 8 

\2 10 11 12 / 



fl 

3 

4 

5 \ 

^5 

1 

3 


A 

9 

6 

7 

00 

^1 

9 

6 

7 

8 


I 2 

10 

11 

12/ 

I 2 

10 

11 

12/ 


/2 1 3 4\ 

4> 5 6 7 8 . (A3) 

\9 10 11 12/ 


To swap the first two particles in the first column, one 
can perform the following sliding operations starting from 
the third-to-last configuration in Eq. (A3): 


/l 3 4 5 \ 

9 6 7 8 

\2 10 11 12 / 


/2 3 4 5\ 

16 7 8 

\9 10 11 12/ 


A 2 3 4\ 

16 7 8 

\9 10 11 12/ 
(A4) 


Equations (A3) and (A4) shows the steps to swap the 
first two particles in the first row, or the first two par¬ 
ticles in the first column, for A = 3 and 5 = 4. This 
can be generalized to any A > 3 and any even 5 > 4. 
The generalization to A > 3 is more obvious because the 
same steps can be directly applied to any A and achieve 
the same goal. The generalization to larger even 5 is less 
obvious. For this case , one needs to repeat the first four 
operations in Eq. (A3) (5/2— 1) times and then perform 
the rest of the steps in Eq. (|A3|) or (|A4|). Since it is possi¬ 


ble to swap any two adjacent particles for any A and any 
even 5 , from the same argument as the A = 5 = 2 case, 
all permutations of particles for any A and any even 5 
are also connected. Similarly, all permutations of parti¬ 
cles for any even A and any 5 are also connected because 
a 90° rotation turns it to the even 5 case. Therefore, all 
permutations are connected as long as A5 is even. 


When A5 is odd, not all permutations are connected. 
This is because none of the sliding operations change the 
parity of the permutation. Thus, two permutations with 
different parity cannot be connected with any combina¬ 
tions of sliding operations. 
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